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Method of predicting shape of wire-like structure , its apparatus 

and its program 

BACKGROUND OF THE INVENTION 

5 The present invention relates to a method of assisting 

a wiring design of a wire-like structure constituted by a 
plurality of line streak members, its apparatus and its program, 
particularly to a method, an apparatus and a program of assisting 
an optimum wiring design of a wire harness wired at a vehicle 
10 as a wire-like structure. 

Normally, at a vehicle or the like, a plurality of 
electrical equipments are mounted thereon, and the electrical 
equipments are connected by a so-to-speak wire harness in which 
15 a plurality of electric wires or communication lines are bundled 
by a bundling member of an insulock or the like or a protective 
member of a tape or the like as the line streak members. As 
shown by Fig. 1, a wire harness 1 is attached with connectors 

2 a, 2 b, 2 c and 2d respective end portions of which are connected 

20 to electrical equipments or the like. Further, various clips 
3a, 3b, 3c, 3d are attached to middle portions thereof , further, 
a branch point 4 is provided thereto. Further, respective 
branch lines of the wire harness 1 constituting from the 
respective end portions to the branch point 4 are basically 
25 provided with different numbers and kinds of constituent line 
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streak members and therefore, the respective branch lines are 
provided with various boldnesses, lengths, elasticities, 
rigidities and the like. Further, in a prior art, a wiring 
design of such a wire harness is frequently carried out mainly 
by hunch and experience of a designer. 

Recently, shape prediction and path design of such a wire 
harness is widely conducted by using a computer support system. 
As a basic technique thereof, a CAD (Computer Aided Design) 
and a general-purpose analysis software such as CAE (Computer 
AidedEngineering) are combinedandutilized. That is, adetailed 
shape of wire harness to which boldness, length type and the 
like are reflected is designed by the CAD, and thereafter the 
shape prediction is calculated by inputting required data by 
number into predetermined general-purpose CAE. After 
evaluating the calculating result. The CADdesigns the predicted 
shape again. Such a cycle is repeated by designers who are skilled 
in operation of CAD and general-purpose CAE. 

[ 0001 ] 

A literature cited in the specification will be shown 
here as follows. 

"Matrix finite element method" by B. Nas, Brain Publishing 
Co., August 10, 1978, p.7-15 

However, the following problem is posed in designing a 
path by repeatedly using CAD and general purpose CAE as in the 



ar t • That i-S/ CAD and general purpose CAE are separated 
in software and therefore, numerical values necessary for 
general purpose CAE need to input at a number of times . Further, 
operation of inputting the necessary numerical values is also 
complicated and enormous time is required for leaning a technique 
of mastering general purpose CAE. Further, according to 
general purpose CAE, there also poses a problem that lead time 
for simulating deformation of a wire harness, that is, a time 
P® ® ^ required forpredicting a shape thereof is mu ch prolonged . 

SUMMARY OF THE INVENTION 

Therefore, in view of the above— described current state, 
it is a problem of the invention to provide a method of predicting 
a shape of a wire-like structure, its apparatus and its program 
for enabling to predict a desired shape firmly in a short period 
of time while dispensing with learning of a special technology 
and skill. 

Iri order to solve the aforesaid object, the invention 
is characterized by having the following arrangement. 

(1) A method of predicting, by utilizing a finite element 
method, a shape of a wiring structure in which the wiring 
structure constituted by a line streak member covered by an 
exterior member is regarded as an elastic body in which a 
plurality of beam elements a linearity of which is maintained 
coupled with each other, the method comprising: 



a drawing step of drawing the wiring structure having 
an arbitrary shape by a manual operation; 

a shape characteristic extracting step of extracting 
shape characteristic from the drawn wiring structure; 

a constraining condition extracting step of extracting 
constraining condition from the drawn wiring structure; 

a material characteristic calculating step of calculating 

ma terial characteristic of the drawn wiring structure based 
on the drawn wiring structure; and 

a predicted shape calculating step of calculating a 
predicted shape which is brought into a physical equilibrium 
state by applying the shape characteristic, the constraining 
condition and the material characteristic to the finite element 
method . 

(2) The method according to (1), wherein a characteristic 
table representing relationships among type, boldness and 
material characteristic of the exterior member is previously 
stored, and 

wherein the material characteristic calculating step 
calculates the material characteristic of the drawn wiring 
structure by applying type and boldness of the exterior member 
calculated based on the drawn wire harness to the characteristic 


table . 



( 3 ) 


Themethod according to ( 2 ) , wherein in the characteristic 
table, longitudinal modulus of elasticity or transverse modulus 
of elasticity of the exterior member is treated as the material 
characteristic of the exterior member. 

(4) The method according to (1) further comprising: 

cil^rm outputting step of outputting an alarm when it 
is determined that the predicted shape which is brought into 
the physical equilibrium state is impossible to be calculated 
based on the wiring structure drawn according to the manual 
operation. 

(5) An apparatus of predicting, by utilizing a finite element 
method, a shape of a wiring structure in which the wiring 
structure constituted by a line streak member covered by an 
exterior member is regarded as an elastic body in which a 
plurality of beam elements a linearity of which is maintained 
are coupled with each other, the apparatus comprising: 

a drawing unit for drawing the wiring structure having 
an arbitrary shape by a manual operation; 

a shape characteristic extracting unit for extracting 
shape characteristic from the drawn wiring structure; 

a constraining condition extracting unit for extracting 
constraining condition from the drawn wiring structure; 
a material characteristic calculating unit for 
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calculating material characteristic of the drawn wiring 
structure based on the drawn wiring structure; and 

a predicted shape calculating unit for calculating 
predicted shape which is brought into a physical equilibrium 
state by applying the shape characteristic, the constraining 
condition and the material characteristic to the finite element 
method . 

(6) A computer readable recording medium storing a program 
which cause a computer to function as an apparatus of predicting, 
by utilizing a finite element method, a shape of a wiring 
structure m which the wiring structure constituted by a line 
streak member covered by an exterior member is regarded as an 
elastic body which has a circular section and in which a plurality 
of beam elements a linearity of which is maintained are coupled 

with each other, the program causing the computer to function 
as : 

a drawing unit for drawing the wiring structure having 
arbitrary shape by a manual operation; 

a shape characteristic extracting unit for extracting 
shape characteristic from the drawn wiring structure; 

a constraining condition extracting unit for extracting 
constraining condition from the drawn wiring structure; 

a material characteristic calculating unit for 
calculating material characteristic of the drawn wiring 



structure based on the drawn wiring structure; and 

a predicted shape calculating unit for calculating a 
predicted shape which is brought into a physical equilibrium 
state by applying the shape characteristic, the constraining 
condition and the material characteristic to the finite element 
method . 


BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a view showing an outline of a total shape of 
a wire harness constituting an obj ect according to the invention . 

Fig. 2 is a diagram showing a relationship between 
representative supporting members for supporting the wire 
harness and constrained degrees of freedom. 

Fig. 3A is a view showing an outlook of the wire harness. 
Fig. 3B is a view showing a behavior of digitizing the wire 
harness of Fig. 3A and Fig. 3C is a view expressing the wire 
harness of Fig. 3A by a beam element and a node point. 

Fig. 4 is a diagram for explaining a degree of freedom 
of the wire harness expressed by the beam element and the node 
point . 

Fig. 5A is a diagram expressing the wire harness by three 
beam elements and Fig . 5B is a diagram showing a state of coupling 
the three beam elements of Fig. 5A. 

Fig. 6A is a view showing a behavior of measuring a moment 
of inertia and a longitudinal modulus of elasticity and Fig. 



6B is a view showing a behavior of measuring a polar moment 
of inertia and a transverse modulus of elasticity. 

Fig. 7 is a diagram showing an example of a characteristic 

table . 

5 Fig. 8 is a block diagram showing an example of a hardware 

constitution according to the invention. 

Fig. 9 is a flowchart showing a main processing procedure 
according to the invention. 

Fig. 10 is a flowchart showing a subroutine of calculating 
10 a material characteristic in the main processing procedure of 
Fig. 9. 

Figs . 11A through HE are views respectively exemplifying 
states of deforming the wire harness in respective processing 
procedures of Fig. 9. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

An embodiment of the invention will be explained in 
reference to drawings as follows. First, an explanation will 
be given of a total shape of a wire harness as a wiring structure 
20 constituting an object of design and representative supporting 
members in reference to Fig. 1 and Fig. 2. Fig. 1 is a view 
showing an outline of the total shape of the wire harness 
constituting the object of design according to the embodiment 
of the invention. Fig. 2 is a diagram showing a relationship 
25 between the representative supporting members for supporting 


8 



the wire harness and constrained degrees of freedom. 

As described above, the wire harness 1 constituting the 
object of design according to the embodiment is attached with 
four connectors 2a, 2b, 2c, 2d both end portions of which are 
connected with electrical eguipments, not illustrated, 
attached with the various clips 3a, 3b, 3c and 3d at middle 
portions thereof and is provided with the branch point 4 . The 
respective branch lines of the wire harness 1 are basically 
Provided with different numbers and kinds of respective’ 
constituent line streak members and therefore, the respective 
branch lines are provided with different boldnesses, lengths, 
elasticities, rigidities and the like. 

The respective connectors 2a, 2b, 2c and 2d are attachably 
and detachably fixed at predetermined positions in accordance 
with fixed positions and mounting directions of opposed side 
connectors on sides of the electrical equipments to completely 
constrain end portions of the wire harness. Further, the 
respective cl ips 3a , 3b, 3c, 3d completely constrain or rotatably 
constrain predetermined portions of the wire harness at 
predetermined positions of cabinets, stays or the like of the 
electrical equipments. 


Here, an explanation will be given of the clips. In the 



clips, basically, there are long hole clips and round hole clips . 
The round hole clip is referred to also as rotational clip and 
is constituted by a base seat portion holding the wire harness 
and a support leg inserted into an attaching hole in a shape 
of a round hole provided at the stay or the like. The round 
hole clip is rotatable around Z-axis (in a direction orthogonal 
to an attaching portion) . 

Meanwhile, the long hole clip is referred to also as fixed 
flip and is constituted by a base seat portion holding the wire 
harness and a support leg inserted into an attaching hole in 
a shape of a long hole provided at the stay or the like. A 
sectional shape of the support leg is constituted by a shape 
of a long hole substantially similar to that of the attaching 
hole. The long hole clip is not rotatable around Z-axis. 

Further, in the long hole clip and the round hole clip, 
there are a corrugated long hole clip and a corrugated round 
hole clip rotatable around X-axes (in a longitudinal direction 
of the wire harness) . Constrained degrees of freedom of the 
respective clips in respective axial directions and around 
respective axes are as shown by Fig. 2. 

In Fig. 2, X-axis, Y-axis and Z-axis correspond to three 
axes orthogonal to each other in a right hand local coordinates 



system at respective node points (or also referred to as nodes) 
the wire harness. Although for example, Z— axis is made to 
coincide with a clip axis, a method of determining the axes 
can pertinently be changed in accordance with a function used. 
Further, in the drawing, constrained degrees of freedom of the 
branch point are also shown for reference. Further, a node 
point on the wire harness arbitrarily set other than the 
above-described constrained points is basically completely 
free although not illustrated here. The constrained degrees 
of freedom are respectively set to the respective nodes prior 
to calculating a predicted path, a reaction force or the like 
as described later. 

Next, an explanation will be given of outlines of assuming 
conditions constituting a premise, a utilized theory and basic 
equations in reference to Fig. 3 through Fig. 6. Fig. 3(A) 
is a view showing an outlook of the wire harness. Fig. 3(B) 
is a view showing a state of digitizing the wire harness, and 
Fig. 3(C) is a view showing the wire harness of Fig. 3(A) by 
beam elements and node points . Fig . 4 is a view for explaining 
degrees of freedom of the wire harness shown by the beam elements 
the node points . Fig . 5 (A) is a view showing the wire harness 
by three beam elements and Fig. 5(B) is a view showing a state 
of coupling three of the beam elements of Fig. 5(A) 


First, according to the embodiment, in utilizing a finite 



element method in designing the wire harness, it is assumed 
as follows. 

(1) The wire harness is assumed to be an elastic body. 

(2) The wire harness is assumed to be coupled with beam 
5 elements. 

(3) It is assumed that linearity is maintained at each 
beam element. 

(4) Asection of the wireharness is assumed to be circular. 
By assuming in this way in this embodiment, application 

of the finite element method to the wire harness which has not 
been carried out, can be carried out. 

According to the embodiment, first, the wire harness is 
digitized. That is, as shown by Fig. 3(A), the wire harness 
1 m which a plurality of electric wires 11 are bundled by a 
protective member of a tape 12 or the like can be regarded as 
a continuous body. Next, as shown by Fig . 3(B), the wire harness 
1 is divided (digitized) into a number of beam elements Cl, 

C2 ' C3 ' -• That is, since the wire harness is like a single 
piece of rope and therefore, the wire harness can be regarded 

to be connected with a finite number of pieces of the beam 
elements . 

Therefore, as shown by Fig. 3(C), the wire harness can 
be represented by coupling the plurality of beam elements Cl, 

C2, C3, - by a plurality of nodes Nl, N2, N3, -. Characteristic 



values necessary for the beam element are as follows, 
length 1 (refer to Fig. 3(B)) 
sectional area A (refer to Fig. 3(B)) 
moment of inertia I 
polar moment of inertia J 
density p 

longitudinal modulus of elasticity E 
transverse modulus of elasticity G 

Further, in the specification, shape characteristics are 
constituted by the length 1 and the sectional area A and material 
characteristic are constituted by the moment of inertia I, the 
polar moment of inertia J, the density p, the longitudinal 
modulus of elasticity E and the transverse modulus of elasticity 
although a description thereof will be given later. 

Further, as shown by Fig. 4, each beam element C (Cl, 
C2 , c3 ' •") is provided with two of a node point a and a node 
point (3. In a three-dimensional space, the node point a is 
provided with 3 of translation components and 3 of rotation 
components and therefore, the node point a is provided with 
a total of 6 degrees of freedom. Further, the same goes with 

the node point p. Therefore, the beam element C is provided 
with 12 degrees of freedom. 


Further, in the drawing. 
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F X j.: force in xi-axis direction of i-th element 

F yi : force in yi-axis direction of i-th element 

F zi : force in zi-axis direction of i-th element 

M x i : moment around xi-axis of i-th element 

M yi : moment around yi-axis of i-th element 

M zi : moment around zi-axis of i-th element 

U X j.: displacement in xi-axis direction of i-th element 

U yi : displacement in yi-axis direction of i-th element 

U Z i : displacement in zi-axis direction of i-th element 

0xi : angular displacement in xi-axis direction of i-th element 

0 yi : angular displacement in yi-axis direction of i-th element 

0 Z i : angular displacement in zi-axis direction of i-th element 

notation a designates the node point on the left side and notation 

P designates the node point on the right side. 
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If the theory of structure for a structure involving a 
large change in shape such as wiring structure, balance equation 
of the finite element method is generally represented as follows : 


2 0 (K + K g )U = P 

where : 

K: rigid matrix 

K g : geometrical rigid matrix 
U: displacement vector 
25 P: load vector 
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However, since the above equation is algebraically 
non-linear simultaneous equation, the equation cannot be solved 
in actual numerical analysis. Therefore, incremental methods 
5 in which a load value is subdivided and subdivisions are 
subsequently added are employed (also employed in case of 
forcibly displacement) . Thus, the balance equation is 
represented by a following incremental form: 

10 (K + Kg)AU = AP - R 

where : 

AU: value of displacement incremental 
AP: incremental load at incremental step 
15 R: compensation vector of incremental vector 

At each incremental section, the balance equation is 
regarded as linear equation and calculated. At this time, 
generated non-balance force (R at the equation) is decreased 
20 to an allowable range by an iteration method in advance to the 
next step. As series of these algorism, well-known method such 
as Newton-Raphson method, arc-length method and the like is 
utilized. 

Incidentally, with respect to a subject of designating 
25 forcible deformation in the case of shape prediction, for example. 
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since omitting the geometrical rigid matrix of the second term 
ori the left side of the balance equation causes a good result 
in many cases, the present embodiment omits it. 

Further, the rigid matrix the first term on the left side 
is a matrix in which elements of the rigid matrix, values of 
coordinates of each of which is rewritten and changed every 
moment according to each incremental step, are converted into 
global coordinates system and are collected. Specific 
representation of element rigid matrix as the basis of this 
is as following equation 1. 


[Equation 1] 









An explanation will be given here of a compatibility 
condition and a condition of equilibrium. Here, for simplicity, 
as shown by Fig. 5(A) , the wire harness is represented by 3 
of the beam elements Cl, C2, C3. In this case, displacements 
of the node point If} of the beam elements Cl and the node point 
2ot of the beam element C2 are equal to each other and forces 
^Pplish to the two node points are also balanced.. Similarly, 

also displacements of the node point 2(3 of the beam element 
C2 and the node point 3ot of the beam element C3 are equal to 
each other and forces applied to the two node points are also 
balanced. Therefore, the beam elements Cl and C2 and the beam 
elements C2 and C3 can be coupled as shown by Fig. 5(B) by 
satisfying continuity of the displacements and the condition 
of equilibrium. 

Further, in the drawing, 

Fxi : force in xi-axis direction of i-th element 

Fyi : force in yi— axis direction of i — th element 

Fzi : force in zi-axis direction of i-th element 

M xi : moment around xi-axis of i-th element 

M yi : moment around yi-axis of i-th element 

M Z i : moment around zi-axis of i-th element 

Uxi • displacement in xi — axis direction of i — th element 



Uyi : displacement in yi-axis direction of i-th element 
U zi : displacement in zi-axis direction of i-th element 
0 xi : angular displacement in xi-axis direction of i-th element 
Gyi : angular displacement in yi-axis direction of i-th element 
0 Z i : angular displacement in zi — axis direction of i — th element 
and i=la, 10, 2a, 20, 3a, 30. 

Further, when the continuity of the displacements and 
equilibrium of forces of the beam elements Cl, C2, C3 shown 
in Fig. 5(B) are shown by a style similar to that of Equation 
(2) as mentioned above. Equation (3) is derived as follows. 


[Equation 3] 



Ml 



2 0 








Here, matrices of 12 rows and 12 columns Ml, M2 and M3 
in Equation (3) are similar to those in shown in Equation (2) . 

5 However, at the portions M12 and M13 at which the matrices Ml, 
M2 and M3 overlap each other, respective constituent elements 
of the respective matrices are added together. 

Further, 4 or more of the beam elements can similarly 
10 be dealt with. In this way, an equation model of a wire harness 
divided into an arbitrary number of beam elements can be formed . 

Incidentally, Equation (3) is simply expressed as 
follows . 

15 [K] {x} = {F} - (4) 

Therefore, a predicted shape can be calculated by 
obtaining respective elements of the displacement vector {x} 
20 which are unknowns based on Equation (3) or Equation (4) . In 
calculating a solution of Equation (3) or Equation (4), a 
well-known arc length method or positive solution method can 
be utilized. Further, the above-described general matrix 
finite element method is introduced in, for example, nonpatent 
25 literature 1, mentioned above. 
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Next, an example of a way of calculating a longitudinal 
modulus of elasticity and a transverse modulus of elasticity 
necessary for predicting a shape according to the invention 
w iil he shown below. Fig. 6(A) is a view showing a behavior 
of measuring a moment of inertia and a longitudinal modulus 
of elasticity and Fig. 6(B) is a view showing a behavior of 
measuring a polar moment of inertia and a transverse modulus 
of elasticity.. Fig. 7 is a diagram showing an example of a 
characteristic table. 

Further, the longitudinal . modulus of elasticity E can 
be represented by Equation (5) , shown below, when a measuring 
method shown in Fig. 6(A) is carried out. 

E=FL 3 / 3X1 - (5) 

Further, the moment of inertia I can be represented by 
Equation (6), shown below, since the wire harness is assumed 
to be constituted by a circular section as mentioned above. 

I=7tD 4 /64 ... (6) 

Therefore, the following equation is established. 

E=64 FL 3 / 3X7lD 4 - (7) 

According to the measurement, a relationship between F 
and x may be measured by constituting the following equation. 

E= (F/X) x (64L 3 /37tD 4 ) 



Meanwhile, according to the invention, the longitudinal 
modulus of elasticity E is calculated for each kind and each 
boldness of an exterior member. That is, as shown by Fig. 7, 
the longitudinal modulus of elasticity E is calculated as 
5 described above for each of boldnesses of an exterior member 
of W/H diameter 1-5, W/H diameter 6-10, W/H diameter 11-15 of 
each kind of exteriors of no exterior, tape-wound 1, tape-wound 

2, — . In Fig. 7, values of 10, 12, ••-, and the like are values 
of the longitudinal modulus of elasticity E of the wire harness 
10 in correspondence with predetermined kind of exterior and-W/H 
diameter . 

In Fig. 7, no exterior indicates a wire harness which 
is not provided with the exterior member of cable tie, tape 
15 or the like. Further, tape-wound 1 indicates a wire harness, 
for example, roughly wound with a tape of a predetermined kind, 
tape-wound 2 indicates a wire harness, for example, half-lapped 
with a predetermined kind of tape. Other than these, various 
kinds of exteriors are assumed by combinations of kinds of ways 
20 of winding of double half-lap, scroll-winding and the like, 
kinds of tapes, wound thicknesses and so on. 

Further, in Fig. 7, W/H diameter 1-5, W/H diameter 6-10, 
W/H diameter 11-15 express boldnesses of respective wire 
25 harnesses by a unit of millimeter. The longitudinal moduli 
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of elasticity E are calculated as described above respectively 
for all of combinations of kinds of exteriors and W/H diameters 
which are realistically feasible to form a data base. 

Meanwhile, the transverse modulus of elasticity E can 
be expressed by Equation (8), shown below, when a measuring 
method shown in Fig. 6(B) is carried out. 

G= ( TL/0 J) x2 - (8) 

The polar moment of inertia J can be expressed by Equation 
(9), shown below, since the wire harness is assumed to have 
a circular section. 

J=7iD 4 /32 - (9) 

Further, torsional force can be expressed as follows. 

T=TS .- (10) 

Therefore, the following equation is established. 

G=(32FSL/0tcD 4 )x2=(F/0) ( 32SL/tiD 4 ) x2 - (11) 

Therefore, the polar modulus of elasticity G can be calculated 
by measuring a relationship between F and 0. 

Each value of the transverse modulus of elasticity G can 
also be calculated for each kind and each boldness of the exterior 
member similar to the longitudinal modulus of elasticity E. 

In this case, each value in Fig. 7 is the transverse modulus 
of elasticity G of the wire harness in correspondence with a 
predetermined exterior kind and W/H diameter. Further, the 



transverse moduli of elasticity G are calculated as described 
above respectively for all of combinations of exterior kinds 
and W/H diameters which are realistically feasible to form a 
data base. 

Incidentally, there is a relationship of Equation (12) , 
shown below, between the transverse modulus of elasticity G 
and the longitudinal modulus of elasticity E. 

G=E/2 ( 1+p) - (12) ..... 

where notation designates Poisson's ratio. 

Therefore, the transverse modulus of elasticity G can 
also be derived from the longitudinal modulus of elasticity 
E and Equation (12) by giving a pertinent Poisson's ratio p. 
Conversely, the longitudinal modulus of elasticity E can also 
be derived from the transverse modulus of elasticity G . Further, 
the above-described measuring methods are only examples and 
respective values of the transverse modulus of elasticity G 
^^d the longitudinal modulus of elasticity E may be acquired 
by methods other than the above -described measurement examples . 


Next, an explanation will be given of a hardware 
constitution according to the embodiment for calculating and 
outputting the shape of the wire harness in accordance with 
a processing procedure, mentioned later, by utilizing the 



above-described theory and basic equations . Fig. 8 is a block 
diagram showing the hardware constitution according to the 
embodiment . 

As shown by Fig. 8, according to the embodiment, there 
is used, for example, a personal computer constituted by 
including a microcomputer 21, an input apparatus 22, a display 
apparatus 23, a printing apparatus 24, a storing apparatus 25, 
a communication interface 26. .and a read/write apparatus 27. 
The microcomputer 21 includes CPU 21a (central processing unit ) , 
ROM 21b for storing boot programs and the like, and RAM 21c 
for temporarily storing various processing results . The input 
apparatus 22 is a keyboard, a mouse or the like for inputting 
the above-described various values and the like, the display 
apparatus 23 is LCD, CRT or the like for displaying the processing 
results and the printing apparatus 24 is a printer for printing 
the processing results. 

Further, the storing apparatus 25 is a hard disk drive 
for storing an installed wiring design assisting program 29a 
and a processing result by the program 29a, the communicating 
interface 26 is a modem board or the like for carrying out data 
communication with an external apparatus by using, for example, 
the internet, an LAN circuit or the like. The read/write 
apparatus 27 is an apparatus of reading the wiring design 



assisting program 29a stored in a record medium 29 of CD-ROM, 
DVD-ROM or the like and writing a calculation result by the 
w ^ r ^- n< =f design assisting program 29a to the record medium 29. 
The respective constituent elements are connected via an inner 
bus 28. 

The microcomputer 21 installs the wiring design assisting 
program 29a read by the read/write apparatus 27 to the storing 
apparatus 25. Further, when a power source is inputted, the 
microcomputer 21. is started in accordance with boot programs 
stored to ROM 21b to start the installed wiring design assisting 
program 29a. Further, the microcomputer 21 carries out 
processing with regard to wiring design assistance according 
to the invention in accordance with the wiring design assisting 
program 29a, outputs a processing result from the display 
apparatus 23 or the printing apparatus 24 and holds the 
processing result at the storing apparatus 25 or the record 
medium 29. The wiring design assisting program 29a can be 
installed also to other personal computer or the like having 
the above-described basic constitution and makes the computer 
function as wiring design assisting apparatus after having been 
installed. Further, the wiring design assisting program 29a 
may be provided not only from the record medium 29 but by way 
of a communication network of the internet, LAN or the like. 


Further, an explanation will be 


given of a processing 



procedure according to the embodiment in reference to Fig. 9 
Fig . 11. Fig. 9 is a flowchart showing a main processing 
procedure according to the invention which is carried out by 
using the hardware constitution shown in Fig. 8. Fig. 10 is 
5 a flowchart showing a subroutine of calculating the material 

characteristic in the main processing procedure of Fig . 9. Fig. 

11(A) through Fig. 11(D) are respectively views exemplifying 
states of deforming the wire harness in the respective processing 
procedures of Fig. 9. . .. 

10 

First, at step SI shown in Fig. 9, a wire harness having 
an arbitrary shape is drawn in correspondence with manual 
operation. The step SI is similar to a designing and drawing 
method using well-known CAD. That is, with respect to a shape 
the wire harness in a straight state on the display apparatus 
23 as shown by la of Fig. 11A, a mouse or the like as the input 
apparatus 22 is used, for example, an end portion thereof is 
displaced from la5 to lb5 as shown by ml. Further, a shape 
of the wire harness having desired length and boldness as 
20 indicated by lb of the drawing is drawn on the premise of being 
arranged in a predetermined three-dimensional space at inside 
of a vehicle (actually with an image as shown by, for example. 
Fig. 3(A) ) . At this occasion, the wire harness is drawn such 
that also the exterior kind explained in reference to Fig. 7 
25 can be discriminated (actually, for example, with an image as 
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shown by Fig. 3 (A) ) . Further, a supporting member of a clip 
or the like is also drawn . Further, in this case, for simplicity, 
it is assumed that only an end portion designated by lal is 
completely constrained and the end portion designated by la5 
5 is deformable . Other portion is assumed to be completely free . 

Further, portions designated by lal through la5 in Fig. 
11 (A) are portions in correspondence with representative node 
points allocated when the wire harness is coupled with a 
10 plurality of beam elements to form a model. Portions ' 
respectively in correspondence with respective portions 
designated by lal through la5 are designated as lbl through 
lb5, lcl through lcS, ldl through ld5, lei through le5 in Fig. 
11(A) through Fig. 11(E). Further, step SI corresponds to 
15 drawing step and drawing unit in claims. 

Successively, at step S2, the shape characteristic is 
extracted from the shape lb of the wire harness drawn at step 
SI. In details, the kind of the exterior member, the boldness 
20 and the length of the wire harness and so on are extracted from 
the shape lb. Further, these are temporarily stored to a work 
area of RAM 21c for processings to be carried out thereafter. 
Step S2 corresponds to shape characteristic extracting step 
and shape characteristic extracting unit in claims. 

25 
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Successively, at step S3, the material characteristic 
is calculated from the shape lb of the wire harness drawn at 
step SI . In this case the material characteristic is 
constituted by the sectional area A, the moment of inertia I, 
the polar moment of inertia J, the longitudinal modulus of 
elasticity E and the transverse modulus of elasticity G included 
in respective elements of the rigidity matrix [K] . Step S3 
corresponds to material characteristic calculating step and 
material characteristic calculating unit in claims. 

In details, in calculating the material characteristic 
at step S3, as shown by Fig . 10, atstepS31, first, the sectional 
area A is calculated. The sectional area A can be calculated 

based on the boldness of the wire harness extracted at step 
S2 . 


Next, at step S32, the moment of inertia I is calculated. 
The moment of inertia I can be calculated in accordance with 
the sectional area A calculated at step S31 and Equation (6) , 
mentioned above. Next, atstepS33, the polar moment of inertia 
J is calculated. The polar moment of inertia J can be calculated 
in accordance with the sectional area A calculated at step S31 
and Equation (9), mentioned above. 


Next, at step S34 , the longitudinal modulus of elasticity 



E is acquired in reference to the characteristic table 25a 
previously stored in the storing apparatus 25. That is, the 
longitudinal modulus of elasticity E in accordance with the 
kind of the exterior member and the boldness of the wire harness 
extracted at step S2 can be acquired easily and firmly. 

Further, the transverse modulus of elasticity G is 
acquired in reference to the characteristic table 25a at step 
S35.. That is, the longitudinal modulus of elasticity E in 
accordance with the kind of the exterior member and the boldness 
of the wire harness extracted at step S2 can be acquired. 
Naturally, as described above, the transverse modulus of 
elasticity G may be derived from the longitudinal modulus of 
elasticity E by utilizing Equation (12) 

When calculation of the material characteristic has been 
finished, successively, at step S4, the constraining condition 
is extracted. The constraining condition is constituted by 
kinds of constraint (completely constraint, rotatably 
constraint, completely free and the like) as shown by Fig. 2 
at the allocated respective node points, coordinates of the 
end portion lbl which is completely constrained and the end 
portion lb5 in the destination of deformation by ml which are 
particularly important. Further, the respective values set 
here relate to the respective elements in the deformation vector 



{x} of the equation (3) . Step S4 corresponds to constraining 
condition extracting step and constraining condition 
extracting unit in claims. 

Successively, at step S5, a predicated shape is calculated 
by using a finite element method. That is, since the shape 
lb of the wire harness shown in Fig. 11 (A) is drawn simply by 
using CAD, the physical equilibrium state is not necessarily 
established to satisfy the shape characteristic, 'the material 
characteristic and the constraining condition and therefore, 
a predicted shape which is brought into the equilibrium state 
is calculated as shown by lc of Fig. 11(B). Step S5 corresponds 
to predicted shape calculating step and predicted shape 
calculating unit in claims. 

In calculation by using the finite element method, as 
shown by Equation (3) shown above, the wire harness is formed 
into a model of coupling a plurality of beam elements and unknowns 
in themodel are calculatedby applying the shape characteristic, 
the material characteristic and the constraining condition to 
the model. That is, when unknowns in the displacement vector 
{x} are calculated, the predicted shape of the wire harness 
in the equilibrium state as shown by lc is provided. 


However, there can be also a case in which a solution 
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of satisfying the shape characteristic, the material 
characteristic and the constraining condition is not provided. 
That is, there can also be a case in which the predicted shape 
brought into the physical equilibrium state as shownby lc cannot 
5 be calculated. By assuming such a case, at step S6, 

success /unsuccess of the predicted shape calculation is 
determined. 

At step S6, when it is determined to have succeeded in 
10 the predicted shape calculation (Y of step S6) , the operation 
proceeds to step S7 and the predicted shape designated by lc 
is drawn in place of the shape designated by lb which is being 
drawn. That is, when the end portion of the wire harness is 
displaced at the input apparatus 22, a shape physically in 
15 equilibrium is automatically calculated by the finite element 
method and the shape is drawn again. That is, the shape 
physically in equilibrium is automatically calculated and drawn 
at the background by following shape design by well-known 
operation by CAD. 

20 

Meanwhile, at step S6, when the predicted shape 
calculation is determined to be unsuccessful (N of step S6) , 
the operation proceeds to step S8 and an alarm stating that 
the predicted shape calculation is unsuccessful is outputted. 

25 The alarm may be displayed on the display apparatus 23 by 
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characters or alarm sound may be emitted by sound emitting 
apparatus (not illustrated) . Further, in this case, the shape 
of the wire harness which is being drawn designated by lb stays 
to be displayed on the display apparatus 23. In this way, when 
the predicted shape physically in equilibrium cannot be 
calculated, the statement is alarmed and therefore, precise 
path design can be carried out in a shorter period of time. 
Step S8 corresponds to alarm outputting step and the alarm 
outputting -unit in claims. - 

Successively, at step S9, presence or absence of a change 
in the shape in correspondence with manual operation is 
determined. For example, assume that the end portion is 
displaced from lc5 to ld5 as shownbym2 of Fig. 11 (c) by operating 
the input apparatus 22 since an obstacle of an electric equipment 
or the like is present on the predicted shape designated by 
lc. Then, the operation returns to step S4 (Y of step S9) to 
extract a new constraining condition. That is, coordinates 
of the end portion lc5 are updated to coordinates of the end 
portion ld5. Further, the operation proceeds to step S5 to 
carry out calculation using the finite element method again 
to thereby calculate a shape satisfying the new constraining 
condition as shown by le of Fig. 11 (D) 


Further, as shown by Fig . 11(E), a shape of the wire harness 


designated by le is drawn on the display apparatus 23 in place 
of the shape designated by Id. Naturally, the processings of 
step S4 through step S8 are repeated when there are changes 
in the shape in correspondence with manual operation again and 
5 again for some reason. When it is determined that there is 
not the change m the shape in correspondence with the manual 
operation (Nof ste P S9) , a series of the processings are finished . 
In finishing, the predicted shape le may be preserved in a hard 
disk or the like as the storing apparatus 25 or may be printed 
10 on paper from the printing apparatus 24. 

Although according to the embodiment, node points other 
than the end portion aremade to be completely free, the predicted 
shape may be calculated in consideration of influence of 
15 gravitational force. The influence of gravitational force 
relates to respective elements in the force vector {F} in 
Equation (3) . 

In this way, according to the embodiment of the invention, 

20 the wire harness is regarded to be the elastic body having the 
circular shape and coupled with the plurality of beam elements 
maintaining the linearity to form the model and therefore, an 
amount of calculation is considerably reduced and the shape 
can firmly be predicted in a short period of time. Further, 

25 data necessary for calculation of the finite element method 
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is automatically calculated from the wire harness drawn in 
correspondence with manual operation and therefore, it is not 
necessary to input numerical values for calculating the 

predicted shape . Therefore, learning of special technique and 
5 skill are dispensed with. 

Further, the method and the apparatus of the invention 
are not limited to the wire harness wired at inside of the vehicle 
but applicable similarly to a wire harness wired indoors. 

10 Further, the method and the apparatus is similarly applicable 
also to a wire harness having a plurality of branch lines branched 
from a trunk line. 

According to the invention, the wire-like structure is 
15 regarded to be the elastic body having the circular section 
and coupled with the plurality of beam elements maintaining 
the linearity to form the model and therefore, in comparison 
with a general finite element model allocated with meshes, an 
amount of calculation is remarkably reduced. Therefore, the 
20 shape can be predicted firmly in a short period of time . Further, 

data necessary for calculation of the finite element method 
is automatically calculated from the wire-like structure which 
is drawn in correspondence with manual operation and therefore, 
it is not necessary to input numerical values for calculating 
25 the predicted shape. Therefore, learning of a special 
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technique and skill are dispensed with. 

According to the invention, the material characteristic 
necessary for calculationof the finite elementmethodcan easily 
5 and firmly be provided from the wire-like structure which is 
drawn in correspondence with manual operation by previously 
referring to the characteristic table in correspondence with 
the relationship among the kind, the boldness and the material 
characteristic of- the exterior- member'. - . .. . 

10 

According to the invention, the longitudinal modulus of 
elasticity and the transverse modulus of elasticity necessary 
for calculation of the finite element method can easily and 
firmly be provided from the wire-like structure which is drawn 
15 in correspondence with manual operation by referring to the 
characteristic table. 

According to the invention, when the predicted shape 
physically in equilibrium cannot be calculated from the 
20 wire-like structure formed in correspondence with manual 
operation, the statement is alarmed and therefore, precise path 
design can be carried out in a shorter period of time. 
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